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START-UP  OF  A  PULSED  BEAM  FREE  ELECTRON  LASER  (FEL)  OSCILLATOR 

A  number  of  successful  free  electron  laser  oscillator  experiments  have 

been  reported. Simple  considerations  concerning  the  spontaneous  radiation 

3  4 

level  Indicated  start-up  times  much  shorter  than  those  observed.  * 

Therefore,  since  a  number  of  experiments  utilizing  shorter  electron  beam 
macropulses  are  being  constructed  or  planned,  there  Is  concern  that  these 
forthcoming  experiments  may  not  be  able  to  achieve  saturation.  A  quantitative 
understanding  of  the  basic  process  governing  the  growth  of  coherent  stimulated 
radiation  from  incoherent  spontaneous  emission  Is  therefore  highly 
desirable.  This  is  especially  true  since  the  device  is  very  sensitive  to 
small  changes  In  parameters.  A  quantitative  understanding  of  the  coherent 
gain  Is  available,  but  the  analysis  of  the  coupling  to  the  incoherent  emission 
is  incomplete.  Here  we  outline  a  classical  one-dimensional  theory  of  the 
start-up  of  the  FEL  oscillator  in  the  cold,  small  signal  regime.  Certain 
important  three  dimensional  effects  are  incorporated  heuristically  by  means  of 
filling  factors.  The  statistical  features  of  the  problem  lead  to  a 
formulation  in  terms  of  an  ensemble  averaged  energy  density  matrix  e,  a 
diagonal  element  of  which  is  proportional  to  the  fraction  of  the 
electromagnetic  field  energy  in  the  associated  Fourier  component.  This  matrix 
obeys  a  linear  equation  in  which  the  inhomogeneous  term  represents  the 
ensemble  average  emission,  and  gain  and  loss  appear  in  coefficient  matrices. 
The  non-diagonal  terms  of  e  yield  information  on  the  cross  correlations 
between  Fourier  components. 

Theories  of  the  FEL^”^  proceed  from  a  continuum  description  of  the 
electron  dynamics,  either  fluid  equations  or  the  Vlasov  equation.  A  proper 
description  of  the  start-up  of  an  FEL  oscillator,  however,  must  take  into 
account  the  fact  the  electrons  are  discrete  and  substantially  uncorrelated, 

since  it  is  the  acceleration  radiation  of  individual  electrons  in  the  wiggler 
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that  provides  the  initial  fields.  The  acceleration  radiation  is  then 
aaplified  by  the  collective  gain  mechanism  associated  with  the  continuum 
description.  Thus  a  statistical  theory  is  required,  couched  in  terms  of 
objects  bilinear  in  the  fluctuating  quantities  so  that  ensemble  averages  are 
non-zero,  even  when  the  ensemble  averages  of  the  fluctuating  current  density 
vanish. 

The  theory  developed  here  is  one-dimensional  in  space  and  treats  the 
electrons  as  governed  by  Maxwell's  equations.  The  wiggler  is  approximated  by 

A 

the  vector  potential  A  *  A  cos(k  z)  0 (z  -  L  )  0 (L  +  L  -  z)  e  where 

0(x)  is  the  step  function  and  Aw  and  !<„  -  27t/1w  are  constants.  The  radiation 

vector  potential  is  written  as 

OO 

,t)  *  £  aR(t)  8*n  NiZ  exP^u>nt)  ex+  c*c*>  where  kn=  wn/c  =  nrr/L,  the 

n  *  1 

separation  between  the  mirrors  is  L  (see  Fig  1),  an(t)  is  taken  to  be  slowly 

varying  in  time,  and  the  tangential  component  of  the  electric 

field  E  *  -  c_13A  /3t  vanishes  on  the  mirrors.  In  what  follows  we  will  assume 
~  ~R 

that  1^1  »  |Ag|.  The  current  density  driving  Ar  can  be  written 

as  J(z,t)  ■  J  (z,t)  +  J.  (z,t)  where  J  is  the  coherent  current  driving  the 
—c  '“lnc  -“c 

stimulated  radiation  (gain)  and  J.  is  the  incoherent  current  due  to  the 

~inc 

discrete  nature  of  the  electrons  and  is  responsible  for  the  spontaneous 

radiation  (shot  noise).  The  coherent  and  incoherent  current  densities  are 

respectively  given  by  J  «  -  |e|v  F  <n(z,t)>  and 

c 

J.  •  -  lelvF  [n(z,t)  -  <n(z,t)>]  where  v  -  eg  -  I e | A  (z)/y  m  c  is  the 
~mc  ~w  me  ~w  ~w  ~w  o  o 

wiggler  velocity.  The  actual  discrete  electron  density  is  n(z,t) 
m  a\>  I  z  “  *j(0)  where  is  the  transverse  electron  beam  area  and  the 
lateral  distribution  of  electrons  has  been  treated  as  uniform.  Only  axial 
discreteness  is  Included.  The  ensemble  average  over  the  initially 
uncorrelated  electrons  is  denoted  by  <  >,  hence  <n(z,t)>  is  the  continuum 
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electron  density  and  n(z,t)  -  <n(z,t)>  represents  the  fluctuating  part  of  the 

density.  The  usual  filling  factor  Is  Fc  -  <?b/or,  where  is  the  transverse 

area  of  the  resonator  radiation  mode;  Flnc  is  the  filling  factor  associated 

with  the  Incoherent  radiation  and  is  somewhat  more  involved.  It  can  be  shown 
—2  -1  —2 

that  F,  -  a. (X,y)  where  X  -  £  (l  +  v  /c)  y  is  the  characteristic 
me  t>  L  o  l  w  oz  qz 

laser  wavelength  and  fm  is  a  loss  factor  due  to  the  finite  size  of  the  mirror 

located  at  z  *  L.  The  loss  factor  for  the  incoherent  radiation  is 

2 

f  *  [2rr  /(I  +  vs  )L1  where  r  is  the  mirror  radius.  In  computing  f  we 
m  1  m  '  w  J  m  m 

have  taken  the  incoherent  radiation  divergence  angle  to  be  **  (1/y  +  8  ). 

w 

The  equation  of  motion  for  the  jth  electron  is 
*j  *  -  ( |e|/Yomoc)2(3/3z  -  c  2v<jz3/3t)(Aw(z)  •  A^z.t))  |  ^  where  the 

right  hand  side  is  the  ponderomotive  acceleration.  The  initial  conditions  are 
*j<tj)  "  0.  ^d  ZjCtj)  -  vqz. 

The  radiation  vector  potential  obeys 
(32/3z2  -  c"232/3t2  -  c"2  vd/dtlA*-  -4vc-1J,  where  v  -  u  /Q, 
u>*  »  2nc/X,  is  the  characteristic  upshifted  frequency,  v  -  c8  is  the 

unperturbed  axial  pulse  velocity,  and  Q  is  the  quality  factor  of  the 
resonator. 

These  equations  when  linearized  yield  a  set  of  coupled  linear  equations, 

the  homogeneous  part  of  which  describes  the  coherent  phenomena,  and  the 

inhomogeneous  part  (involving  the  fluctuation  in  the  electron  density) 

represents  the  incoherent  emission.  The  incoherent  (spontaneous)  emission 

part  of  a.  satisfies  <a  >  ■  O.Thus  in  order  to  obtain  non-trivial  statistical 
n  n 

information  one  defines  the  total  radiation  energy  density  matrix 

e  (t)  ■  k  k  <a  (t)  a  (t)>.  Note  that  e  -  ’  2  <|a  |2>  is  proportional  to 
nm  n  m  n  m  nn  n  n 

the  ensemble  average  energy  density  in  the  nth  Fourier  component,  and  e  is 

nm 

related  to  the  correlation  of  the  nth  and  mth  Fourier  components  of  the 
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electric  field.  These  statistical  quantities  are  measurable.  It  may  be  shown 
after  considerable  calculation  that  when  the  gain  per  pass  is  small 


e(tN+  t)  -  (1  -  —■)  e(tN)  +  G  (tN,  T)  •  e(tN)  +  e(tN)  •  GH(tN,i)  +  S(tN,T) 
c 

(1) 


where  the  matrix  G  enbodies  the  coherent  response,  S  represents  the 
spontaneous  emission,  t^  is  the  time  the  Nth  electron  pulse  entered  the 
wlggler  and  H  denotes  the  Hermitian  conjugate.  Indeed  the  trace  of  S  is 
proportional  to  the  total  spontaneous  energy  radiated  by  one  pulse  of 
electrons  in  traversing  the  wiggler.  The  elements  of  the  gain  matrix  and  the 
source  matrix  at  time  tN  +  x  are  given  by 


(tN«) 


0  a.  a. 

It  _b  2iir(n-m)(N-l)  5L/L  „ 

32  L  Y  2  V  c 

o  c 

2  2  2 

*b  *  novw  Fjnc  2ix(n-m)(N-l)6L/L 
2L  Lo^  6  anm 


nmpnm8nm^T  ^ 


P  h  (t) 
nm  nm 


» 


where 


3  ix  x  i(x-x) 

8nm(T) - 2^e  ”(l  +rL)  8lnxn  '  e 

xx  n 

n  m 


sin(x  -  x  ) 
( - 


x  -  X 
n  m 


x 

ra 


2ix 

e 


h  (t) 
nm 


TT  ( ( sin2x„+  sin2xm  -  sin2(xn  -  xj) 


-  i(sinx  cosx  -  sinx  cosx  -  sin(x  -  x  )  cos(x  -  x  )H, 
nn  mm  nm  'n  m  11  * 
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*„«  •  V<lH,o)CL.)*  V  'W  ckn<1  '  vo J^n- 


nm 


✓f  -UVW4) 

— r—  e  ,  for  Gaussian  electron  pulse  profile 


sin  (k  -  k  U. /2 
n  o  D 

(k  -  k  )  lj'2  ~ 
n  m  ft 


,  for  square  shaped  electron  pulse  profile. 
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0  «  v  /c,  w.  ■  4n|e|  n  /m  and  6L  ■  L  -  L./20  .  The  matrices  e  and  S  are 

o  o  "b  o  o  bo  a  a 

Hermit lan. 

Equation  (1)  has  been  integrated  numerically  to  obtain  c.  The  parameters 
employed  in  the  numerical  studies  are  given  in  Table  1,  corresponding  to  the 
Stanford  experiment.^  Figure  2  shows  the  peak  radiation  power  within  the 
resonator  as  a  function  of  the  number  of  beam  pulses  that  have  traversed  the 
resonator  for  various  values  of  resonator  length  mismatch  6L  «  L  -  L,  /20  . 

D  O 

Figure  3  shows  the  asymptotic  gain  as  a  function  of  6L.  The  mirror 

_3 

mismatch  61  *  -  1.1  x  10  cm  corresponds  to  maximum  gain  but  not  maximum 

saturated  power.  Maximum  saturated  power  occurs  for  6L  between  0  and  -l.l  x 

—3  -3 

10  cm.  The  range  in  6L  for  nonzero  gain  is  -  3.0  x  10  cm<  6L  <  0,  in  fair 

_3 

agreement  with  the  experimental  range  of  2.5  x  10  cm.  The  maximum 
calculated  multi-mode  (finite  beam  pulse)  energy  gain  is  0.16  whereas  the 
single  mode  (continuous  beam)  yields  a  value  of  0.25.  Finite  beam  pulse 
effects  therefore  reduce  the  linear  gain  by  approximately  60Z.  The  maximum 
observed  gain  is  0.10. 

Figure  A  shows  the  spatial  distribution  of  the  electron  pulses  (square) 

and  the  radiation  power  pulse  at  the  entrance  and  exit  of  the  wiggler 
-3 

for  6L  ■  -  1.0  x  10  cm.  Upon  entering  the  wiggler  the  radiation  pulse 
slightly  lags  the  beam  pulse,  while  at  the  exit  of  the  wiggler  the  two  are 
completely  overlapped.  Figure  5  shows  the  energy  spectrum  of  the  radiation 
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pulse  after  100  beam  pulses*  The  spontaneous  radiation  energy  spectrum  Is 

also  shown  for  reference  purposes. 

Equation  (4)  suggests  that  one  can  roughly  compute  the  relationship 

between  P„,  the  peak  power  in  the  resonator  after  the  Nth  pulse,  to  P  ,  the 
N  o 

power  emitted  spontaneously,  by  assuming  a  constant  average  gain  per  pass  g. 

An  elementary  calculation  yields  when  N  »  1  and  g  «  1,  Pn/Pq  “  N  -  i  + 

N 

(1  +  g)  ■  N  +  exp  (gN).  Clearly  when  gN  »  1  the  result  is  very  sensitive  to 

small  changes  in  g  and  N.  If  one  takes  the  experimental  values  corresponding 

to  the  maximum  observed  final  power  of  *  2.7  x  10^  W  within  the  resonator, 

-2 

N  «  540  and  the  computed  spontaneous  power  of  Pq  •  6.5  x  10  W,  one  finds 
that  g  *  0.037.  The  experimental  value  of  the  linear  gain  is  0.067.  In  view 
of  the  sensitivity  to  changes  in  N  and  g  the  results  are  not  inconsistent. 
Moreover  this  effective  value  of  g  is  smaller  than  the  linear  gain  predicted 
by  the  present  model  which  is  reasonable  since  non-linear  effects  and  initial 
beam  thermal  effects  must  lower  the  gain.  Unfortunately  the  currently 
available  data  is  inadequate  to  make  other  detailed  comparisons  with  this 
small-signal  theory. 

Our  analysis  immediately  suggest  possible  ways  to  substantially  shorten 
the  oscillator  start-up  time  while  maintaining  high  saturated  power  levels. 

The  first  approach  takes  advantage  of  the  fact  that  the  maximum  linear  gain 
and  maximum  saturated  power  occur  for  different  values  of  6L,  which  we  will 
respectively  denote  by  and  By  slightly  increasing  the  frequency  of 

the  R.F.  accelerating  field,  “>acc>  during  the  start-up  period,  i.e., 
decreasing  the  beam  pulse  separation,  the  value  of  5L,  could  be  varied  from  an 
initial  value  of  6L^  to  the  value  of  SL^,  thus,  decreasing  the  start-up  time 
while  maintaining  high  final  power  levels.  The  required  fractional  increase 
In  uacc  is  |5L^-  10~®  for  the  parameters  of  ref.  (3,4).  The  same 
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effect  may  also  be  realized  by  simply  changing  (increasing)  the  mirror 
separation  during  the  start-up  periods  Another  possible  method  of  decreasing 
the  start-up  time  would  be  to  simply  increase  that  part  of  ^inc  associated 
with  mirror  losses,  i.e..  Increase  fm«  This  could  be  accomplished  by 
increasing  the  effective  size  of  the  mirror  located  at  z  ■  L.  The  additional 
extension  of  the  mirror  would  necessarily  have  a  different  curvature.  This 
last  approach  should  make  it  possible  to  contain  a  far  larger  portion  of  the 
incoherent  radiation. 
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Beam  Energy,  (y  -  l)m  c 
o  0 

43  MeV 

Total  Gamma ,  yq 

85 

Axial  Gamma,  Y 

02 

69 

Peak  Current ,  1 

P 

1.3  A 

Pulse  Width, 

0.75  mm 

Pulse  Separation,  L 

D 

25.4  m 

Beam  Radius,  r. 

D 

0.25  mm 

Wlggler 

Parameters 

Wavelength, 

3.3  cm 

Amplitude  (helical),  B^ 

2.3  kG 

Length,  Lw 

5.3  m 

Resonator 

and  Radiation 

Resonator  Length,  L 

12.7  m 

Resonator  Losses  (round  trip) 

1.5% 

Radiation  Wavelength,  X^ 

3.3  pm 

Spot  Size,  r 

o 

0.167  cm 

Beam  Filling  Factor,  F 

c 

0.017 

Incoh.  Rad.  Loss  Factor,  fD 

0.05 

RESONATOR 


Fig.  1  Schematic  of  a  pulsed  electron  beam  FEL  oscillator. 
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Fig.  2  Peak  power  of  the  filiation  pulse  as  a  function  of  the  number  of 
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Fig. 


RADIATION  POWER  (RELATIVE  UNITS) 


Fig.  4  Radiation  pulse  power  relative  to  the  spatial  distribution  of  the 

electron  pulse  (square)  at  the  entrance  of  the  wiggler  (t  -  t„)  and 

N 

exit  of  wiggler  (t  -  tN+  Lw/vzo>»  where  N  »  1  denotes  the  electron 
pulse  number. 
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ENERGY  (RELATIVE  UNITS) 
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